Signatures of integrability in charge and thermal transport in ID quantum systems 
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Integrable and non-integrable systems have very different transport properties. In this work, we 
highlight these differences for specific one dimensional models of interacting lattice fermions using 
numerical exact diagonalization. We calculate the finite temperature adiabatic stiffness (or Drude 
weight) and isothermal stiffness (or "Meissner" stiffness) in electrical and thermal transport and 
also compute the complete momentum and frequency dependent dynamical conductivities a(q,ui) 
and K(q,u>). The Meissner stiffness goes to zero rapidly with system size for both integrable and 
non-integrable systems. The Drude weight shows signs of diffusion in the non-integrable system and 
ballistic behavior in the integrable system. The dynamical conductivities are also consistent with 
ballistic and diffusive behavior in the integrable and non-integrable systems respectively. 

PACS numbers: 72.10 -d, 72.10 Be 



A perfect metal is distinguished from a superconduc- 
tor by the absence of the Meissner effect, although both 
share an infinite conductivity at zero frequency, fn one 
dimension, another fundamental kind of metallic systems 
is possible, corresponding to an integrable system. We 
term this as an integrable metal (IM) and such a sys- 
tem can display characteristics that are different from 
both a perfect metal and a superconductor, such as an 
infinite conductivity at all temperatures. In contrast, a 
perfect but non integrable metal (NIM) becomes resis- 
tive at finite temperatures due to inelastic collisions and 
umklapp. We track this difference between an IM and 
NIM directly in this work by calculating various trans- 
port stiffnesses also the dynamical conductivities for both 
charge and energy using numerical exact diagonalization 
and a suitably generalized Kubo formula The main 
result is that the NIM shows evidence of diffusive behav- 
ior in both the static and dynamic conductivities, while 
the IM shows several signs of ballistic transport. 

The adiabatic stiffness (or Drude weight) is given by [1( 
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Here T a = — lim^o r [J a (k), K a {— k)], where K a (k) is 
the Fourier component of a local density operator (a — e 
for the charge density, a = E for the energy density 
and a = Q for the heat density) while J a (k) is that of 



the corresponding current operator. Here p n 
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is the Boltzmann weight of the many-body state with 
energy e„ and the partition function. This is arises in 
the adiabatic evolution of an equilibrium state disturbed 
in the infinitely remote past by the application of a time- 
dependent perturbation, and is the coefficient of the delta 
function in the conductivity for charge transport. 

Another stiffness, the isothermal (or "Meissner") stiff- 
ness D a can be defined by the same expression as Eqn. Q] 
by omitting the Kronecker delta function that forbids 



equal energy states. For charge transport, this arises in a 
study of the Byers-Yang type curvature of the free energy 
with respect to a flux through the ring On general 
grounds, a superconductor which displays the Meissner 
effect has a non-zero value of D e . While a non- vanishing 
D a implies a non- vanishing D a , the converse is not true. 

The distinction between D a and D a resides in the rel- 
ative statistical weight of states with equal energy. It 
also follows that if D a is known to vanish, then D a can 
be written as a sum over only equal energy states 
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The current operator commutes with the crystal momen- 
tum, so n and m have the same value of crystal momen- 
tum q. In general, an IM, can have degeneracies among 
states with same q due to the presence of dynamical sym- 
metries. We can thus write D as the sum of D^, the con- 
tribution from states degenerate with others at the same 
value of q and -D™ d , the contribution from states with no 
degeneracies within the same crystal momentum sector. 

In an interesting paper it has been argued that in 
an IM with time reversal invariance, a non-zero value of 
D e is due to the presence of a large number of current 
carrying states, which are degenerate with a partner state 
at the opposite value of q. An important question is how 
these current carrying states contribute to and D™ 1 ; 
What is the role of the dynamical symmetries in deter- 
mining the magnitude of the adiabatic stiffness? This is 
especially relevant when time reversal is broken by e.g. 
an irrational flux through the ring, as used here, whereby 
even non-degenerate states can carry currents. 

We have been motivated by these questions to under- 
take an exhaustive numerical study of a typical IM and 
NIM, by computing their exact energy spectra and all the 
current matrix elements and consequently the transport 
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stiffnesses and dynamic conductivities. This exercise is 
done in ID with a popular model, the t — t' — V model 
of spinless fermions, where a single non-zero parameter 
t' destroys the integrability, but lets the system remain 
a perfect metal (an NIM). All stiffnesses are generally 
non-zero for finite systems, and it is only the systematics 
of their size dependence that gives us reliable informa- 
tion on the large size behavior. We also study for the 
first time, the complete dynamic conductivity a(q, u>) and 
thermal conductivity n(q, u>) in an IM and NIM, using an 
appropriate scheme for binning the discrete data points 
to obtain continuous functions. 





FIG. 1: (Top) The isothermal and adiabatic (inset) stiffness 
for electrical transport in the integrable system at various val- 
ues of (5 with an irrational flux twist. It can be seen that the 
former tends to zero rapidly with increasing system size while 
the latter appears to go to a constant, indicative of ballistic 
transport. (Bottom) The same stiffnesses for electrical trans- 
port in the non-integrable system again with an irrational flux 
twist. Both tend to zero with system size but D e does so more 
rapidly. 

The t — t' — V model of spinless fermions on a ID ring 
is given by the Hamiltonian 



+c]c ]+2 ) + Vj2( n j - l/2)K+i - 1/2). (3) 



For concreteness, we choose the value of t and V to 
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FIG. 2: Plots of a(q, oj) (above) and n(q, uj) (below) vs. ui for 
the integrable system with L — 16 at f3 = 0.001. It can be 
seen in the plot of n(q, oj) that there is a band of frequencies 
for a given value of q in which the value of n(q, ui) is large and 
then abruptly falls to zero at the boundaries. This behavior 
is much less pronounced in a(q,u>). 



be 1.0 and 2.0 respectively. With t' = 0, this model 
can be mapped onto the integrable spin 1/2 Heisenberg 
ring using the Jordan- Wigner transformation. We set 
t' = 1 when we model the NIM. Chains of length 6-18 
are studied and we focus on half-filling for definiteness. 
The Hamiltonian Eqn. [3] is then numerically diagonal- 
ized exactly. We would like to emphasize that our con- 
clusions vis a vis integrability and non-integrability are 
quite general and independent of the microscopic param- 
eters and the filling. We also apply an irrational flux twist 
to break time-reversal invariance, which too does not af- 
fect the integrability. For convenience, we will hence- 
forth use the term non- degenerate only for states which 
are not degenerate with any other state at the same q. 
With the irrational flux twist, the NIM now has only 
non-degenerate states while the IM has both degenerate 
and non-degenerate states. This has the advantage of en- 
abling us to analyze the behavior of the IM in contrast to 
the NIM, mainly in terms of degeneracies (and set Df = 
for the NIM) . We calculate D a and D a for electrical and 
thermal transport at different temperatures (T = 1/(3 is 
measured in units of i). 

We show here only the results for electrical transport 
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FIG. 3: Colored contour plots of a(q, u) (above) and K(q,u>) 
(below) for the IM with L = 16 at f3 = 0.001. Momentum 
(q = to 8). An eighth order polynomial interpolation has 
been used to convert the discrete numerical data into a func- 
tion of a continuous momentum variable q and the feature 
marked A on the lower plot is an artifact of that. The plot 
of n(q, uj) shows the banded continuum in the center bounded 
by dark blue regions. 
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FIG. 4: Plots of a(q, uj) (above) and n(q, uj) (below) vs. uj for 
the NIM with L = 16 at /3 = 0.001. n(q, uj) no longer displays 
the banded behavior of the integrable system and falls off 
gradually (and not abruptly) at higher frequencies. The peaks 
in both a(q,ui) and K,(q,u) at large q are also present in the 
density-density correlation functions indicating the presence 
of overdamped excitations. 



in Fig. [T] but the ones for thermal transport display ex- 
actly the same qualitative behavior. The issue of whether 
the charge stiffness (D e in our model) in the Heisenberg 
ring (t' — 0) is actually non-zero at finite temperature 
remains contentious 0,0- Our results up to the largest 
system size (L = 18) seem to indicate that it is non-zero. 
Indeed, the overall picture that emerges is consistent with 
ballistic transport. A confirmation of this to rule out any 
slow (logarithmic) decay with L is currently beyond the 
range of accessible system sizes. The NIM on the other 
hand shows a clear (exponential at high temperature) 
decay of D e with system size indicating diffusion. This 
behavior persists to temperatures higher than the typi- 
cal energy scales in the problem and in fact even up to 
infinite temperature. These results are interesting since 
they show that the system is aware of its integrability 
or non-integrability even up to infinite temperature. D e 
vanishes rapidly in the thermodynamic limit in both sys- 
tems (even without the flux twist). This shows in partic- 
ular that the IM despite having an infinite conductivity 
at uj — is not a superconductor. It is a curious fact 



that D e is positive at small system sizes in the IM and 
negative in the NIM. The situation is reversed without a 
flux twist. The transport stiffnesses in small systems can 
in general have either signs depending on the microscop- 
ies p. 

We now outline a statistical study of D e in the two 
systems. Since the Meissner stiffness D e is zero at large 
L, D e is given by Eqn. O To enable a comparison of 
D e in the IM and the NIM at a given L, we remove the 
dependence on microscopic parameters (i.e. normalize) 
by dividing by T e . We also normalize Df and D™ d in 
the same way. For the largest system sizes (L = 18), we 
find that the normalized D e is about 8-15 times larger 
in the IM than the NIM for the values of (3 considered. 
One can now ask the question we set out to answer in 
a slightly different from: Are the large current carrying 
states which give a large (normalized) D e in the IM com- 
pared to the NIM in the degenerate states in the former, 
or in the non-degenerate ones? This can be answered 
by comparing the magnitudes of the normalized and 
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Motivated by this question, we have studied the statis- 
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tics of the number of degeneracies and the the current 
matrix elements in the two systems. We find that at 
L = 18 in the IM, about 10% of the states are pairwise 
degenerate with a very negligible fraction of higher order 
degeneracies. Further, most of the degeneracies occur 
among states with total lattice momentum or ir. The 
ratio of the nrmalized Df/D™ d is at most 0.2 for the val- 
ues of (3 investigated. On the other hand, the ratio of 
£>™ d between the IM and the NIM is about 8-15. Thus, 
we come to the following interesting conclusion: Even 
though the IM has degeneracies, their contribution to 
the charge stiffness is not significant. The large current 
carrying states are primarily among the non- degenerate 
states. We also note that the eigenvalues of the cur- 
rent operator in a given degenerate sector in the IM, 
are always roughly of the same magnitude. We reiter- 
ate that the term non-degenerate is used here only for 
states with no degeneracies within a sector of total mo- 
mentum q. Thus, it appears the presence of dynamical 
symmetries does not have a direct effect on the charge 
stiffness through the creation of degenerate states. A 
more detailed account of the statistics alluded to above 
will be presented elsewhere. 

We now focus on transport at finite frequency (uj) and 
q. The q and uj dependent conductivities are given by 

A a (q,uj) = c(u>) ^ Pn\(n\J a (q)\m}\ 2 5(e m ~e n -huj), 

(4) 



where c(uj) = j- 
A E (q,uj) = n(q,uj) 



A e (q,uj) = a(q,uj) and 



These conductivities can also be 
related to density-density correlation functions of the 
charge and energy [f|. We choose a small value of 
(3 = 0.001 since for a calculation of this sort, numerical 
exact diagonalization is most efficient only at very high 
temperatures [7(. The irrational flux twist turns out to 
be unimportant to the results here. o~(q, uj) and n(q, uj) as 
functions of uj (in units of the hopping t) for different val- 
ues of q (going from to 8) for an L — 16 IM are shown in 
Fig. H n(q = 0,uj) cx S(uj) since [J E (q = 0),H] = but 
[J e (q = 0),H] ^ and thus a(q — 0,w) has some struc- 
ture at uj 0. A more interesting feature is that n(q, uj) is 
non-zero only within a band of frequencies for small val- 
ues of q and goes to zero abruptly at the boundaries of 
the band. Moreover, this band shifts to higher frequen- 
cies with increasing q. The dispersion of these banded 
modes appears to be roughly linear (at small q), consis- 
tent with considerations of ballistic transport. The situa- 
tion is similar to the case of free fermions where ballistic 
transport causes a similar banded structure in n(q,uj). 
The common feature of this IM and free fermions is inte- 
grability, which it appears is strongly associated with the 
concept of ballistic transport. o~(q, uj) does not display 



the same banded feature as n(q, uj) as prominently, pre- 
sumably due to the fact that [J e {q = Q),H] ^0, where 
the analogy with free fermions does not apply. Fig. [3] 
shows contour plots of a(q, uj) and n(q, uj) to better illus- 
trate the banded nature of n(q, uj) and enable comparison 
to cr(q,uj). 

For the sake of comparison we also present numerical 
data for a(q, uj) and K,(q, uj) of the NIM. The plots of these 
quantities are shown in Fig. |U [Jb(q = 0), H] ^ here 
and n(q — 0,uj) is non-zero at uj ^ 0. Both o~(q = 0,uj) 
and n(q — 0,uj) display a finite non-analytic singularity 
at uj = 0, which has been attributed to diffusion and non- 
linear hydrodynamics for the former In this system 
n{q,uj) does not display the banded behavior of the IM 
and gradually goes to zero with increasing uj for all q. 
This is indicative of diffusion in this system, which we 
have verified by also computing the density-density cor- 
relators directly. The features (bumps) at large values of 
q(= 5 — 8) and uj also appear in these correlators indicat- 
ing the presence of interesting overdamped excitations, 
which will be investigated in detail elsewhere. 

To conclude, we have demonstrated that the IM shows 
several signs of ballistic behavior as opposed to diffusion 
in the NIM. This has been illustrated through an exten- 
sive numerical calculation of finite temperature trans- 
port stiffnesses and dynamical conductivities. Futher, 
the contribution of degenerate states coupled by the cur- 
rent has been shown to be insignificant to the Drudc 
weight in the IM. The NIM seems to possess overdamped 
excitations at large momenta 
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